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Abstract 

We propose a simple stochastic model describing the position-dependent ball-passing net¬ 
work in football games. In this network, a player on a certain area in the divided fields is 
a node, and a pass between two nodes corresponds to an edge. Our model is characterized 
by the consecutive choice of a node dependent on its intrinsic fitness. We derive the explicit 
expression of the degree distribution, and find that the derived distribution reproduces the 
real data quit well. 
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1 Introduction 


In the past years, scientific studies on football have attracted growing interest jlj. It has been 
suggested that there are some statistical laws in football dynamics, including goal distributions 
MSI, temporal features of the ball touches self-similarity of the movement of the ball and 
players [6]. 

Complex network analysis mi is another approach to extract statistical properties from 
football games. In particular, a network in ball passing is composed by nodes and edges cor¬ 
responding to players and passes, respectively. Several works have proposed the assessment 
methods of players and teams based on some network measures such as clustering coefficient, 
betweenness centrality, and PageRank mm- Furthermore, the structural properties and the 
spatiotemporal patterns of ball-passing networks have been also investigated mm- 

Previously, we have proposed a method to create a “position-dependent” ball-passing network 

in football games [13] - In this network, a node represents a player on a certain area in the divided 
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fields, and an edge corresponds to a pass between two nodes as shown in Fig. [T] Note that a 
node is distinguished by the combination of who makes or receives a pass in which area; that is, 
one player defines different nodes according to the area of a pass. The degree k in this network 
corresponds to the sum of the numbers of making and receiving passes of each node. We have 
obtained the degree distributions of this network from real data. It was found that they were 
fitted well by a truncated-gamma distribution, whose probability density function is given by 

f(k)=Afk^ 1 e-^ (1) 

where the domain of k is given as 0 < k < k mguX , and u, A, and fc max are the fitting parameters. 
The normalization constant M depends on &: max . By introducing fe max , we obtained the better 
fitting results than the ordinary gamma distribution. 

In our previous paper |13| . we have also proposed the numerical model called the “Markov- 
chain model” describing the ball-passing process. In this model, the passing sequence is simu¬ 
lated by a Markov chain. The ball-possession probability of node j at time t is calculated 
by the Markov chain 

N 

(t+i) _ (t) p . 

i =1 

where Pi^j is the transition probability of a ball from node i to j, and N is the total number of 
the nodes. We have assumed that x® is proportional to the degree of node j. For the football 
games, Pi~>j is assumed to be defined as 

D _ Qi r ij) x P{Lj) 

~ Z 

where Q(rij ) represents the rate of completing a pass dependent on the distance between 

the two nodes i and j. R(Lj ) denotes the factor for existence probability of the player receiving 

a pass, where Lj is the distance of the node j from its home position. (The home position is 

assigned for each player randomly.) Z is the normalization constant which is determined to 

satisfy Pi->j =1- We have discussed numerically the property of the cumulative distribution 
j 

G(x) of the ball-possession probabilities Xj, as a substitute for the degree distribution, and 
compared G(x) with the truncated-gamma distribution Q. We have found that G(x) depends 
mainly on the factor R(L), and reproduces the truncated-gamma distribution by choosing an 
appropriate function for R(L). Meanwhile, in the simplified condition where Pi^j oc R(Lj), we 
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get = R(Lj). Then, the probability distribution g(x) of x^\ defined as g(x ) = —dG(x)/dx, 
and the probability density function p(L) of Lj hold the relation, g(x)dx = p(L)dL. And, g(x) 
can be expressed as 


g(x) = p[R l {Zx)\ — [i? 1 (Zx)] 


(2) 


On the basis of the above results, the present paper derives the explicit expression of the 
degree distribution based on the above Markov-chain model. For the derivation, we employ the 
framework of the fitness model, and extend it to the case where networks contain the temporal 
feature, i.e., the time series of the passes. 


y 




Fig. 1. A schematic representation of the creation of a position-dependent ball-passing network. 
Each player on the three areas in the left figure correspond to the three nodes in the right 
figure, and the passes between these players correspond to the edges. Each node is given the 
serial numbers, such as “159”, “106”, and “131”, which are determined from the area and player 
numbers. Further information about the creation of the network is summarised in our previous 
paper [13J. 


2 Extended Markov-chain model for ball passing 

2.1 Setup 

In this section, we propose an extended model based on the above Markov-chain model. In 
the extension, we take into account the interaction of two teams A and B, each of which consists 
of 11 players. A node expresses a player on a divided area in the field. The field is divided into 
2A sections along the goal direction, and A sections along the direction vertical to the goal 
direction. (The field division shown in Fig. [T] is the case A = 3.) Thus, the total number of 
divided areas is 2A 2 , and each area equally has 22 nodes; nodes are distributed uniformly over 
all areas. The number of nodes in each team denoted by Na and TVg, is 11 x 2A 2 , and the total 
number of nodes are N = Na + Ng. Each node in team A and B is given the serial numbers a 
(a = 1,... Na) and b (b = Na + 1,..., Na + Ng), respectively. For each player, we assign one 
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of areas in the divided field as the home position, and define the distance L between the home 
position and the position for each node. L is the quantity defined for all nodes individually, 
and we call L a “fitness” hereafter. (Fitness is a term usually used in complex network analysis 
mm .) The probability distribution of L for each team is expressed as pa(L) and ps(L). 

Now, the passing sequences are assumed to be random transfer of the ball between nodes. 
For the transition probability between two nodes, we assume the following four forms depending 
on the teams to which these nodes belong: 


Pa'^a — P aRa(B a) / ZAi 

(3a) 

Pb’^b = PbRb^P^IZb, 

(3b) 

Pa’^b = (1 — Pa)Rb(L} } )/Zb , 

(3c) 

Pb'^a = (1 ~ Vb)Ra(Lo)/Za- 

(3d) 


Here r/A and t]b represent the ball-passing probability within the same teams, and the exis¬ 
tence probability R(Lj) of node j is a monotonically decreasing function of L. Owing to the 
normalization of Pi^-j 

N N 

yi Pa'- >j — 1’ Pb’-tj ~ 1) 

3 = 1 3 = 1 

the coefficients Za and Zb are expressed as 

N a N 

Za = ^ Ra(Lo,' ), Zb = ^ Rb{Lv). 

a '—1 b'=N a+1 

The explicit forms of pa{L) and Ra(L) are discussed in Sec. 3. 


Team A 




> 1 A 


Fig. 2. A schematic representation of the extended model. Teams A and B have Na and Nb 
nodes respectively, and each node has the fitness L. 
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2.2 Ball-possession probability 

We introduce the ball-possession probability x ( p (j = 1,, Na, Na +1,..., N) for the node 
j at time t, and the probability vector am' is defined as 

T (t) _ r it) (t) it) it) it) (*)-] 


a; 


^ is normalized at each f as 


N a 

Ed> = W 

a= 1 


TV 


E 4 l) = W 

6=7V a +1 


where xj^ and are the ball-possession probabilities for team A and B which satisfy + 


xg> = 1. 


The time evolution of Xa' 1 is given by the Markov chain 


N 

T (t) -SX r (t-l)p. 
x a ~ / J x ? - r j^m- 

3 = 1 


(4) 


Substituting Eqs. (3a) and (3d) into Eq. Q, we obtain 


(t) (*—!) RA(L a ) n 

= X a' VA Za + 2^ x\, ; (1 - 7)b) 


X\’ = 


a '=1 


b'=JV A +l 


Ra(Lo) 

Za 


(1 - Vb) + ( T)A + VB- l)X\ 


it- 1) 


Ra(Lo) 

Za 


(5) 


By taking summation for a from 1 to Na in the both sides of Eq. (©, the following recurrence 


relation for X A is obtained: 


X A = (4 - vb) + (VA + VB- l)^f 1] - 


The solution is 


do) 


x«‘> = - + ^ + + _ I)t 

2 - rj A ~ VB 2 - rjA - r) B 

Substituting Eq. ([6]) into Eq. ([5]), we have 


T (t) _ 


1 - VB (1 - f?s) + (va + VB - 2)X^ 


( 0 ) 


■( va+vb -1 y 


Ra{L 0 ) 
Za ' 


( 6 ) 


(7) 


2 - VA ~ vb 2 -r)A-VB 

Here, the second term of Eq. 0 decreases exponentially with time t, because \va + Vb ~ M < 1 
holds from 0 < t]a, Vb < 1 - The total number T of ball transitions in one game is more than 
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1000 (see Table [ 2 ] in detail), so the exponential term becomes ignorably small compared with 
the first constant term. In the stationary state, the ball-possession probability for node a in 
team A can be written as 

_ 1 ~ rj B Ra(Lo) 

2 — 7? A -VB Za 

The same discussion is applied to team B. Therefore, we can write the ball-possession probability 
for a node with fitness L in team A and B as 


x a (L) = 
x B (L) = 


1 -VB Ra{L ) 


2 - r] A - t)b Z A 
1 -VA Rb(L ) 


2 -r) A - VB Z B 

2.3 Derivation of the degree distribution 


Here, we focus on a node a in team A whose fitness is L, and derive the in-degree and out- 
degree distributions f^ n \k') and /^ out ■*(£;”). The in-degree of node a increases when the node 
receives a pass from another node in the same team. This event takes place at each time step 
with probability 

N A 

^ X al P a i^ a = VAXa(L). 
a '=1 

On the other hand, the out-degree of node a increases by the event that the node makes a pass 
to another node in the same team, and it occurs with probability 

N A 

^2 x aPa->a> = VAXa(L). 
a'=l 

Here, we introduce the conditional probability q^ CL \k'\L) (qA Ut \k”\L)^ that a node with fitness 
L gets the in-degree k' (out-degree k”) at time T. The in-degree or out-degree of the node with 
fitness L is increased by one at each time step with the probability t)axa{L). Since these events 
are independent each other, q^'' > (k'\L) and q^ ut \k” \L) are given as the binomial distribution: 

L) = (TW.4(£)]*'[ 1 - VA m(£)] t -‘\ 

sSrVlO = 
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where T is the total number of ball passing. If field division A is large enough, L can be regarded 
as the continuous variable, so that we calculate f A n \k') and f A ut \k”) as follows: 

k') = [r] A XA(L)} k '[l - rjAX A (L)} T - k ' p A (L)dL, 

f^ nt \k”) = [r] A x A (L)] k ' [1 - r] A x A (L)] T ~ k ”p A (L)dL. 

Next, we derive the degree distribution f A {k). Since the degree k is the sum of the in-degree 
k' and the out-degree fc”, the probability q A (k\L ) that the node with fitness L has the degree k 
is given by the convolution 


q A (k\L)= £ q f\k'\L) q ^\k”\L) 

k=k'-\-k” 
f e 2T\ 

k )[VAXA(L)} k [l-VAX A (L)} 2T - k . 


Hence, the degree distribution f A (k) is obtained as follows: 

/*oo 

f A {k) = / q A (k\L)p A (L)dL 

Jo 


lo 


2 T 


[tiaxa{L)T{ 1 -t]axa(L)} 2T k pA{L)dL 


f 

Jo 


°o / 2T 

k 


\Ra(L)] 

k 

i - 

Ra(L) ] 

[ Z'a \ 


Z' A \ 


-t 2 T-k 


p A (L)dL, 


where Z’ A is defined as 


_ 2 - rjA ~ VB y 

£ A — 77 ~ ZjA' 


JA Va{ 1 - Vb) 

By the variable transformation u = Ra(L) / Z' A , fA(k) is rewritten as 

/ “max /O r T\ p 

( k ) [u] k [l - u] 2T - k p A [Rp 1 ( Z’ A u )] — [R-/ (Z' A u)] du, 


du 


where 


_ R A ( oo) _ 

— ryl — U - 




R A ( 0) 

^A 


( 8 ) 


(9) 


From the saddle point method, the following approximation of Eq. © is obtained when 2 T 3> k 
(see Appendix for derivation): 


f A (k) ~ -p A [. R A l {Z' A k/2T)\ — [. R A l {Z' A k/2T)\ . 


( 10 ) 
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3 Application to real data 


To compare our theoretical result with the real data, we have to fix two functions Ra(L ) and 
Pa(L) in Eq. ©■ Since each player is less likely to visit areas distant from its home position, 
Ra(L ) corresponding to the existence probability of each player should be a nronotonically 
decreasing function. Here, we assume 

Ra(L) oc exp[— (L//3)"»], (11) 


where /3 is the scale parameter and m is the shape parameter. We adopt this form because it 
can express some different decreasing functions by changing m; it becomes Gaussian function 
for m = 1, and exponential function for m = 2. 

Next, we determine the function p A (L), which corresponds to the number of nodes whose 
distance from its home position is L. Since such a node is on a circle centered at its home 
position with radius L, and nodes are assumed to be distributed uniformly over all areas, pa{L ) 
grows linearly for small L, but turns to decrease for large L because of the finite field size. The 
following function 


p A {L) = -^Lex p — (L/uj) 2 


UJ 


( 12 ) 


satisfies these properties. The factor exp[—(L/ cj) 2 ] expresses the effect of the field boundary. We 
have checked whether this function is good approximation of the probability distribution of L (see 
Fig. I for a typical example). And also, we have found that u becomes almost constant for any 
configurations of the home position as follows. In the condition where A = 25, we assigned one 
of areas to each player randomly as the home position, and examined the probability distribution 
of L for 1000 samples. For each sample, we fitted the probability distribution of L by using Eq. 


(12), and obtained oj = 22.72 ± 1.34. This result implies that almost the same value of oj is 
obtained irrespective of the configuration of the home position. 
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L 


Fig. 3. (i) A typical example of the configuration for the home position of each player in the 
50 X 25 field division, (ii) The probability distribution of L in this home position. The solid 
curve is the fitting result by Eq. (12). 


Finally we derive the explicit expression of the degree distribution for the ball-passing net¬ 
work. Substituting Eqs. © and ( |12[ ) into Eq. and introducing the parameter p = 

(w//3) 2 / m , we get 


fA(k ) = / 

'J Umin 


2 T 
k 


[u] k [l-u] 


2T-k m /log(Uma x /u) 


d \ 




m— 1 


U 


exp 


log ( U mSLX /u) 
d 


where 


du , 
(13) 


= 0 


Z' 


1 r) A {\-r) B ) 
Z a 2- rj A - r] B ' 


Now, we analyze the real data. We have investigated 18 networks obtained from 9 real games 
as shown in Table [lj with the field division A = 3. Degree distribution [Eq. (13)] contains the 
six parameters, /u, m, T, N A , r] A and r]B- The value of N A is 198 when A = 3. T is evaluated by 
directly counting the number of ball transitions in a whole game, and r) A and tjb, the ball-passing 
probability between the same teams, are obtained as the ratio of the number of ball transitions 
within the same teams and T. These values are summarised in Table [2] Then, we regard the 
remaining two parameters /i and m as the control parameters for fittings. Figure [I] shows the 
cumulative degree distributions of each network in a single logarithmic scale. The solid curves in 
each panel are the cumulative distribution F(k ) of Eq. ©. The values of p, and m for fittings 
are summarised in Table [3j We find that all real data in Fig. [4] are in good agreement with Eq. 


( 13 ). 
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Table 1. Real game data for the analysis. 



Game 

Place 

Date 

Score 

Competition 

(i) 

Japan vs Paraguay 

South Africa 

2010.06.19 

0-0 

Wcup 2010 

(ii) 

Japan vs Vietnam 

Japan 

2011.10.07 

1-0 

Kirin Challenge Cup 

(hi) 

Japan vs Tajikistan 

Japan 

2011.10.11 

8-0 

WCup Asian qualifier 

(iv) 

Japan vs North Korea 

North Korea 

2011.11.15 

0-1 

WCup Asian qualifier 

(v) 

Spain vs Italy 

Poland 

2012.06.10 

1-1 

Euro 2012 

(vi) 

Germany vs Holland 

Ukraine 

2012.06.13 

2-1 

Euro 2012 

(vii) 

Mainz vs Hertha 

Germany 

2013.09.28 

1-3 

Bundesliga 13-14 

(viii) 

Manchester City vs Everton 

England 

2013.10.05 

3-1 

Premiere league 13-14 

(ix) 

Manchester United vs Tottenham 

England 

2014.01.01 

1-2 

Premiere league 13-14 


Table 2. The parameters obtained from the real data. 



Game 

T 

VA 

VB 

(i) 

Japan vs Paraguay 

1457 

0.51 

0.58 

(ii) 

Japan vs Vietnam 

1345 

0.65 

0.55 

(hi) 

Japan vs Tajikistan 

1438 

0.69 

0.39 

(iv) 

Japan vs North Korea 

1094 

0.59 

0.55 

(v) 

Spain vs Italy 

1471 

0.72 

0.62 

(vi) 

Germany vs Holland 

1390 

0.65 

0.69 

(vh) 

Mainz vs Hertha 

1138 

0.47 

0.62 

(viii) 

Manchester City vs Everton 

1314 

0.65 

0.60 

(ix) 

Manchester United vs Tottenham 

1215 

0.61 

0.52 
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(iii) (iv) 
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(viii) 
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(ix) 


Fig. 4. Fitting for the cumulative distributions of real data by using Eq. (13). The panels (i)- 
(ix) correspond to the nine games in Table [lj The solid curves show the cumulative distribution 
function of Eq. (THl . 


Table 3. The values of n and m for fittings of each network. 


Game 

Team 

h 

m 

(1) 

Japan 

3.27 

1.26 

Paraguay 

4.62 

2.25 

(ii) 

Japan 

3.92 

1.10 

Vietnam 

4.70 

1.80 

(hi) 

Japan 

4.90 

1.40 

Tajikistan 

4.20 

1.50 

(iv) 

Japan 

3.97 

1.14 

North Korea 

3.40 

1.80 

(v) 

Spain 

4.36 

1.17 

Italy 

4.57 

1.53 

(vi) 

Germany 

4.35 

1.52 

Holland 

4.40 

1.70 

(vii) 

Mainz 

4.67 

1.71 

Hertha 

5.80 

3.00 

(viii) 

Manchester City 

4.00 

1.20 

Everton 

3.89 

1.80 

(ix) 

Manchester United 

4.00 

2.00 

Tottenham 

3.80 

2.40 


4 Discussion 

Regarding Eq. ( |13[ ), a more simplified expression can be derived. Since k/{2T) < 0.01 from 
the data in Tableland Fig. [4j the following approximation of Eq. ( |13[ ) holds: 

,, , Ttl / log (fcmax 

f A {k) ~ — - 1 - 

A* \ M 


, v \ m— 1 

/k)\ i 


exp 


log (. k max /k ) 




(14) 
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(see Appendix for detail discussion). Here, A; max is the expected value of the maximum degree 
given by 

^ = (15) 

Z A 2 — T] A — tjB 

In Fig. [5j we compare the two functions (13) and (14). It is found that the two functions take 
almost the same values except a high-degree part. Such a mismatch in the high-degree part is 


also found in the results of fittings of the real data by Eq. (14). However, if we also regard 
the maximum degree as another fitting parameter represented by femax, all degree distributions 


are fitted well by Eq. (14) as shown in Fig. |6l The values of £ and the maximum degree 


of real data are summarised in Table. [ 4 J Although kinal and fcma? are different 


m some 


teams such as Paraguay in the game (i) and Hertha in the game (vii), we conclude that Eq. (14) 
roughly describes the degree distributions of ball-passing networks. 

Eq. (14) is the Weibull distribution where we set y = log ( kmax/k ): 

m— 1 


fA(y) = -r- 

/i \n 


exp 


Here m and n correspond to the shape and scale parameters, respectively. (Eq. (14) has been 
used for the statistics in earthquake mm for instance.) We note that Eq. becomes the 

i—i 

power-law function /^(fc) ~ when m = 1, which corresponds to A —> 00 in the truncated- 
gamma function [Eq. 01- In the real data, some teams fulfill this condition; Japan in the game 
(i) , (ii) and (iv), Spain in (v) and Manchester City in (viii) have the value of m close to 1 (see 
Table [3]). This is why the truncated-gamma distribution also fits to the real data. 



Fig. 5. Comparison of cumulative degree distribution given by Eq. (13) (the solid curves) and 
Eq. (14) (the points), where Na = 198, ija = Vb = 0.5, y = 4.5, and m = 1.5. 
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(ix) 


Fig. 6. Fitting for the cumulative distributions of real data by using Eq. (14) (the solid curves). 


Table 4. The values of maximum degree. 


Game 

Team 

JAP 

ft max 

, (real) 
^max 

(i) 

Japan 

20.0 

22.0 

Paraguay 

75.0 

42.0 

(ii) 

Japan 

33.9 

32.0 

Vietnam 

36.5 

29.0 

(hi) 

Japan 

72.0 

63.0 

Tajikistan 

15.5 

14.0 

(iv) 

Japan 

24.4 

23.0 

North Korea 

35.0 

21.0 

(v) 

Spain 

45.0 

42.0 

Italy 

37.0 

29.0 

(vi) 

Germany 

44.0 

35.0 

Holland 

47.0 

44.0 

(vii) 

Mainz 

29.1 

21.0 

Hertha 

500.0 

68.0 

(viii) 

Manchester City 

31.1 

29.0 

Everton 

40.5 

30.0 

(ix) 

Manchester United 

45.0 

39.0 

Tottenham 

18.4 

25.0 


In the present study, we have derived the explicit expression of the degree distribution [Eq. 


(13)] from the extended Markov-chain model. In contrast, our previous study |13| heuristically 
adopted the truncated-gamma distribution for fitting. Although both functions fit to the real 


data well, Eq. (13) has the following three advantages over the previous one. First, Eq. (13) is 
derived from the model based on the simple passing process. The key point of our model is that 
Pi^rj depends only on t]a, ijb, and R(Lj). This is a considerable simplification of the real football 
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games. Nevertheless, Fig. [4] shows that such simplification even preserves essential features of 
the actual degree distributions. The second point is the number of control parameters. The 
truncated-gamma distribution [Eq. 0 ] needs to control the maximum degree fc max to reproduce 
the real data. On the other hand, Eq. (|13|) does not contain the control parameter corresponding 


to the maximum degree since the domain of k of Eq. (13) is given as 0 < k < oo. Thus, Eq. (13) 
can reproduce the real degree distributions with the only two parameters ^ and m, which are 
fewer than the truncated-gamma distribution. Third, fj, and m have the clear physical meaning 
as follows. As shown in Sec. 3, cj is calculated as u = 22.72 ± 1.34 in the case where A = 25, so 
that n ~ (22.7Namely, both fi and m are determined from the function Ra(L). We 
illustrate the dependence of Ra{L) on fi and m in Fig. [7J This figure shows that the variance of 
Ra{L ) decreases with /x, and increases with m. Since Ra(L) represents the existence probability 
of a player dependent on the distance L from its home position, /x and m control the typical 
moving range of the player around its home position. In this sense, /x and m reflect the strategy 
or playing style of a team. However, we emphasize that the degree distribution does not change 
greatly by such differences of each game. We expect that Ra{L) is determined directly by the 
analysis based on the more detailed positional data of each player. 



L 



L 


Fig. 7. Change of the shape of Ra{L) with /j, and m. (i) ^-dependence with m = 1.5, (ii) 
??i-dependence with [i = 4.5. 


From the viewpoint of complex network, our model is related to the fitness model mm- 
The difference of the two models is as follows. In the original fitness model, the algorithm 
of the network creation is static, and that an undirected network with un-weighted edges is 
generated. On the other hand, our model contains the temporal features, i.e., the time series 
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of the passes, and generates a directed network with weighted edges. Hence, our model is not 
only the extension of the original fitness model, but also a kind of a temporal network [18]. The 
ball-passing networks in football game and the networks generated from our model seem to be 
similar to the co-occurrence network for human language mm in that the following method 
is in common use; a node represents a word, and an directed edge connects the two adjacent 
words in the same sentence. In creation of sentences, this type of co-occurrence network reflects 
the process of consecutive choice of words and its time series. Moreover, in such a creation 
process, each word is considered to be chosen according to its importance which is associated 
with the fitness. Therefore the co-occurrence network is expected to modelled by our Markov- 
chain model. We believe that our model makes a theoretical contribution to network science, as 
well as the practical analysis of football games. 

5 Conclusion 


In the present paper, we have extended the previous model describing the position-dependent 
ball-passing network in football games. In the extended model, we have taken the effect of 
the opponent team into account, and assumed that the transition probability depends on the 
existence probability R{Lj) of each node. The explicit expression of the degree distribution (13) 
which has two parameters /z and m are derived, and it can reproduce the real degree distributions 
quite well. We have also found that /z and m are determined from the function R(L), and they 
characterize the moving range of a player in each team. Furthermore, we have shown that 


Eq. (13) is simplified to the function (14) when 2 T 3> k is satisfied, and this simple function 
also approximates the real distribution well. Although our model is simple, it incorporates the 
essential features of the formation of ball-passing networks in the actual football games. 
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Appendix: Derivation of Eq. (10) from Eq. 



Eq. ([9]) can be written as 


f A (k) 



( 16 ) 


where 


<p(u) = klnu + (2 T — k ) ln(l — u), 


d 


h(u) = p A [R-/ (Z' A u)} - [. R- a 1 (Z»] 


tp(u) has the peak at u = k/2T. Expanding ip(u) to the second order near u = k/2T, and 


substituting it into Eq. (16), we get 


f A (k) ~ - ( | exp 


<P 


2 T 


exp 


4T 3 

k(2T-k) [ U ~ 2 T 


h(u)du. 


(17) 


When 2 T 3> k, the integrand in Eq. © decreases rapidly without u = k/2T, and it is allowed 
to expand the integration range to (— 00 , 00 ). Hence, we can approximate Eq. © as follows: 


f A (k) ~ - ( | exp 


2 T 


exp 


<P 




k_ 
2 T 

_k_ 
2 T 


)}f 


exp 


4r3 

k(2T-k) [ U ~ 2T 


du h —— 
2 T 


k(2T — k)ir ( k 

4T3 h \2f 


(18) 


Here, we apply Stirling’s formula, t\ — \/2irt t*e t , to Eq. (18), and use 

_ k k (2T — k) 2T ~ k 


exp 




k 

2T 


(2 T) 


2 T 


we obtain Eq. (10): 


f A {k) ~ 


-Hh) 

-PA { R~a 1 (Z' A k/2T)] 4 [R -1 (Z' A k/2T)} . 
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